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The control of nuclear spins in quantum dots is essential to explore their many-body dynamics 
and exploit their prospects for quantum information processing. We present a unique combination of 
dynamic nuclear spin polarization and electric-dipole-induced spin resonance in an electrostatically 
defined double quantum dot (DQD) exposed to the strongly inhomogeneous field of two on-chip 
nanomagnets. Our experiments provide direct and unrivaled access to the nuclear spin polarization 
distribution and allow us to establish and characterize multiple fixed points. Further, we demon¬ 
strate polarization of the DQD environment by nuclear spin diffusion which significantly stabilizes 
the nuclear spins inside the DQD. 

PACS numbers: 72.25.Pn 73.63.-b, 03.67.-a, 73.63.Kv, 


In III-V semiconductors the weak hyperfine interaction 
between nuclear and electron spins has a strong impact 
on the electron spin dynamics owing to the fact that each 
conduction band electron interacts with a large number 
of nuclei mi- This situation can give rise to dynamic nu¬ 
clear spin polarization (DNSP) [SHHl and exciting many- 
body quantum physics such as complex hysteretic dy¬ 
namics [DISlIIS], multistabilities [a IS], collectively en¬ 
hanced transport and dissipative phase transitions 
[M IS). However, the thermal fluctuations of nuclear 
spins, even present at cryogenic temperatures, also cause 
decoherence of spin qubits [201 El]- Ignoring correlations, 
the influence of the nuclei on the electron spin dynamics 
is usually described within a semiclassical mean held ap¬ 
proach, which expresses the nuclear spin polarization in 
terms of an effective magnetic held for the electron spin, 
Bnuc(r), the Overhauser held [22]. In a DQD charged by 
one electron in each dot, thermally huctuating nuclear 
spins result in a held difference AB„uc in the order of 
few mT between the two dots jH]. In equilibrium, the 
time averages of Bnuc(r) and ABnuc vanish. This huc¬ 
tuating held nevertheless causes a weak mixing of singlet 
and triplet states [S] EH El] being explored for quantum 
information processing I1I1E5]. To control this mixing, 
it would be necessary to stabilize the difference of the 
effective magnetic helds in the two dots, ABnuc- 

Here, we combine DNSP with EDSR to study the dy¬ 
namic polarization of nuclear spins on the one hand, and 
the decay of the polarization on the other hand. We 
demonstrate the existence of multiple attractive hxed 
points (FPs) in the steady-state solution of the driven 
system, where the decay of Bnuc(r) is exactly canceled 
by its dynamical build up 011ES]- Our results demon¬ 


strate that the FPs differ from each other by their spatial 
distributions of Bnuc(r)- Relevant for spin qubit applica¬ 
tions, the singlet-triplet mixing of each FP can thereby 
be fine tuned by adjusting external parameters such as 
the external field Hext [direction as in Fig. [T][a)] or the 
energy detuning e between the two dots [i.e. the singlet 
configurations (1,1) and (2,0), where {n,m) denotes the 
number of electrons in the left (n) and right (m) dot]. At¬ 
tractive FPs in DNSP are often characterized by a narrow 
nuclear spin distribution, hence sharply reduced nuclear 
spin fluctuations, which provides possible advantages for 
the preparation of coupled dots for quantum informa¬ 
tion applications [lEZ]. Our measurements allow us to 
characterize FPs by their nuclear spin distribution, and 
their dynamic and static stability. In particular, we show 
that the diffusion of nuclear spins outside the DQD has a 
strong influence on the build-up and decay dynamics of 
their polarization inside the DQD which it further stabi¬ 
lizes. Such an enhanced stability of FPs promises positive 
impact on the coherence of the electron and nuclear spin 
states. 

Our DQD design, presented in Fig. [^a), incorporates 
two single-domain nanomagnets. They generate an addi¬ 
tional, static inhomogeneous field Bnm(r), such that the 
total effective field difference is ABeff = ABnuc + ABnin. 
In equilibrium, the static jABnmj — 45 mT [28] exceeds 
the fluctuations of ABnuc of ~ 2 mT |3] by far, degrading 
ABnuc to a weak perturbation. Important for qubit ap¬ 
plications this stabilizes the singlet-triplet splitting and 
yields advantages in controlling the nuclear spin dynam¬ 
ics [8]. Furthermore, our sizable ABnm causes a cor¬ 
responding separation of the Zeeman energy in the two 
dots and allows us to resolve EDSR experiments in the 
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FIG. 1: (a) Scanning electron microscope image of a DQD 

defined in a two-dimensional electron system 85 nm beneath 
the surface of a GaAs/AlGaAs heterostructure [28]: GaAs 
surface in gray, gold gates in yellow, cobalt gates in blue. Red 
circles indicate approximate quantum dot positions, green ar¬ 
rows the physical current direction for /tr > /tl , a white arrow 
the direction of Bext. (b) Charge stability diagram /(lA, Vr). 
An arrow marks the detuning axis, a dotted line £ = 0. A 
double arrow indicates pulsing of gate voltages lA. and Vr 
during an EDSR experiment [see labels in (a)]; the modula¬ 
tion VV. = Vf2 + vsm{ujt) with v ~ 3 mV is applied at the red 
point labeled EDSR. (c) /(Bext) as Bext is swept. Variations 
suggest the existence of multiple FPs. (d) Distribution of cur¬ 
rent traces /(Bext) corresponding to FPs I-IV and “others” 
(referring to traces which did not fit to types I-IV) detected 
within 384 sweeps. 


individual dots |55|. We operate our DQD in the vicin¬ 
ity of the (1,1) <->■ (2,0) charge transition. In response 
to the applied voltage V = (^r — ^i)/e = ImV elec¬ 
trons tunnel one by one through the DQD via the cy¬ 
cle (1,0)—>-(1,1)—>■(2,0)—>-(1,0). We measure the resulting 
current which is, however, strongly suppressed by Pauli 
spin blockade (PSB) of the transition (1,1)—>•(2,0) ^5] . 
In the stability diagram of our DQD in Fig. Sb), an 
extended region of PSB is clearly visible as reduced cur¬ 
rent at the base of the current-carrying double triangle. 
Details are explained in Ref. for the identical sample. 

We follow two complementary approaches to study 
the nuclear spin dynamics. First, we actively polarize 
the nuclear spins by sweeping Bext and driving electrons 
through the DQD. We measure the background leakage 
current still flowing in PSB, which grows with increasing 
singlet-triplet mixing being proportional to the compo¬ 
nents of ABeff IS] (if we ignore weak influences of co- 
tunneling and spin-orbit interaction). In our second ap¬ 
proach, we let Bnuc(r) (produced as described above) 
decay and directly measure it by performing electric 
dipole induced spin resonance (EDSR). We have experi¬ 
mented with several scenarios, but for better comparabil¬ 


ity here we discuss only measurements taken under the 
following conditions: we start with an initialization time 
^init = 180 s with Bext = 1/ = 0 at £ ~ 0.1 meV (red dot 
in Fig.j^b) labeled ’’Drag”; the interdot tunnel coupling 
is tuned to tc — 20 peV) to let any remaining Bnuc(r) 
decay. Next, we apply V = 1 mV (at otherwise identical 
settings) and sweep Bext at the rate 0.67mT/s to a h- 
nite value and then sweep back to zero at —3.3 mT/s. To 
avoid complications by long time memory effects, before 
each measurement series we preconditioned the system 
with a number of identical sweeps. 

Figure [^c) shows /(Bext) during typical sweeps to 
B™^’' = 0.6 T. Even though we keep dot and sweep 
parameters identical we find four different character¬ 
istic current traces each one occurring multiple times 
[Fig. Qd)] but in arbitrary order. Within each type 
7(Bext) is reproducible, even including sudden current 
steps and the noise level, see Fig. [^c). The different 
curve types /(Bext) are stable over many minutes (jumps 
between them occur very rarely at sizable Bext, not more 
than in one out of 50 sweeps). They strongly depend on 
Bext, and persist even for very slow sweeps, but are lost if 
the sweep is performed too fast. We will show below, that 
the magnetic field sweeps are accompanied by a build-up 
of sizable (and type-dependent) nuclear polarizations. 

Knowing that / scales with ABes, this all points to 
corresponding stable FPs with steady-state nuclear spin 
configurations Bnuc(r,Bext)- Traces II-IV contain sud¬ 
den jumps (within less than a second), which might indi¬ 
cate reproducible transitions between some FPs at spe¬ 
cific values of Bext- Before each sweep we let Bnuc(i’) de¬ 
cay by setting P = 0, hence / = 0, where DNSP is absent 
(and only the trivial FP in equilibrium with Bnuc(r) = 0 
is left). Our experiments indicate that within the first few 
seconds of / 7 ^ 0 the nuclear spins arrange themselves at 
one of the FPs (even at Bext = 0). Figurej^d) presents a 
statistics of the rate of different curve types. It suggests 
that for our settings four FPs are almost equally likely 
occupied under these conditions. Interestingly, the order 
at which the different types occur is random, likely being 
related to random fluctuations. In our case these could 
be thermal fluctuations of B„uc(r, Bext) |1[30] or random 
telegraph noise in the local DQD potential, called charge 
noise [?TH55] . 

Type I curves [see Fig. [^c)] are characterized by a 
reduction of the current from / ~ 1 pA at Bext — 0 to 
/ < 80 fA at Bext > O.IT. In Ref. [HI type I has been 
associated with the resonance between the T+-triplet and 
the singlet state (both for one electron in each dot). In 
our system, the strong reduction of /(Bext) implies that 
especially the components of ABeff perpendicular to Bext 
should be reduced to AB/^ ~ 0. Likewise, the other FPs 
with larger currents correspond to more inhomogeneous 
Bnuc(i') at sizable Bext- 

Utilizing EDSR experiments, we can measure the ef¬ 
fective magnetic field value to monitor Bnuc(i’). Here, 
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FIG. 2: (a) Measuring scheme and an example cnrrent trace 

I{t) (type III); Bnuc(r) begins to decay at t = 0. (See main 
text for details.) (b) I{t) during the decay of Bnuc(r) for vari¬ 
ous EDSR frequencies. Cnrrent peaks marked with red arrows 
indicate EDSR resonances. (For clarity the upper traces are 
vertically offset in steps of 50 fA.) (c) Initial current decay 
for various t„. (d) EDSR resonance times tpeak tpeak S'® 
fnnction of tw. Horizontal lines are mean values. 


we measure its decay in the two dots after we have 
built up Bnuc(r) by DNSP. As sketched in Fig. [^a), 
we start by sweeping B^xt to within tsweep as de¬ 

scribed above. There, we hold Bext and V constant dur¬ 
ing thoid = 120 s. For the remainder of the experiment 
we keep B^xt = Next, at t = 0, we initiate the 

decay of Bnuc(r) by switching V off to stop DNSP. After 
the waiting period tw, we go back to F = I mV. If tw was 
long enough, the FP is lost and Bnuc(r) continues decay¬ 
ing. To measure this, we perform EDSR by periodically 
pulsing to e ~ — 0.5meV [see Fig. 0b)] and applying 
an rf-modulation at a hxed frequency to the gate voltage 
V^. Details of our EDSR procedure are explained in Ref. 
HSJ At £ <C —tc the electrons, confined in the DQD and 
affected by EDSR, can be considered localized in the in¬ 
dividual dots. As long as the rf-bursts are off-resonance, 
I remains small, but whenever hf = g^BB^s in one of the 
two dots, PSB is lifted and we expect enhanced current. 

Typical traces I{t) are plotted in Fig. 0b) for the 
case of type III curves for three different frequencies and 
= 0.65 T. The frequency independent gradual de¬ 
crease of / within the first 250 s is related to the decay 
of B„uc(r) away from the FP. It indicates that the FP is 
close to a singlet-triplet resonance, where I has a max¬ 
imum. On top of this background, we observe current 
spikes indicating EDSR resonances [arrows in Fig.j^b)]. 
Interestingly, we find two distinct resonances per decay 
curve, the first occurring at and the second at 
which we assign to the right and left dot, respectively (see 
below). A resonant EDSR experiment can also lead to 



FIG. 3; (a) EDSR resonance frequency versus time /(tpeak): 

ipeak (blue, right dot) and tpeak (green, left dot) for type I (III) 
curves in the upper (lower) panel. Lines are fits according to 
Eq. 0. (b) Type III decay curves /(tpeak) for various 
(c) Slowing down the decay of the nuclear spin polarization: 

tpeak (thold ). 


DNSP (more so at higher frequencies) [HI fTHl IHlUHi)] . The 
uppermost I{t) curve in Fig. 0b) measured at / = 2 GHz 
demonstrates this effect, where the second EDSR spike 
repeats multiple times and causes / to fluctuate strongly. 
Important for our analysis, no DNSP is induced as long 
as / is off-resonant. This is always the case in each dot 
before the corresponding current peak occurs the first 
time, marking the relevant tpeak and ^pg^i^. 

For measuring the undisturbed decay of Bnpc(r), it 
would be desirable to avoid any DNSP effects during the 
decay. As a continuous F = 0 and no rf-modulation are 
no option, we carefully monitor and set up the experi¬ 
ments to avoid unwanted DNSP effects. In particular, 
we found that the background current does not cause 
any DNSP away from the FPs: for instance, variations 
of tw (at F = 0) do not influence the initial decay of 
/(t), see Fig. 0c). In addition, and are in¬ 
dependent of w as demonstrated in Fig. 0d). However, 
a correlation is evident in the jitter between tpgai. and 
tpeak (the correlation coefficient is 0.87). Such a behav¬ 
ior is expected if the origin of the fluctuations is related 
to noise occurring with frequencies small compared to 
l/(^peak ~ ^peak)' ^his favors charge noise, which has 
its weight at long time scales |3IM1I37], over the ther¬ 
mal fluctuations of Bnuc(r), which are fast compared to 
^peak ^peak m- Charge noise modulates the geometry 
of the confinement potential. The details of the latter 
determine parameters with a possibly strong impact on 
the static and dynamic properties of the FPs, including 
the detuning e, the positions of the charge centers, con¬ 
sequently ABjiin, and the decay of Bnuc(r). 

To map out the decay of Bnuc(r) in the two dots, we 
present tpgak tpgak (x-axis) of type I and HI curves 
for various frequencies (y-axis) and = 0.36 T in 
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Fig. [^a). The solid lines are theory curves, where we as¬ 
sume that Bnuc(r) is antiparallel to Bgxt (an assumption 
fully compatible with our data 0) and decays exponen¬ 
tially, such that 


hf = \g^^B Kr + Km + Kuc(t = 0) exp(-t/T)]| , (1) 


where we replaced vectors by their z-components (along 
the external held), a good approximation for the data in 
Fig.^a ,b) |33]. In accordance with EDSR measurements 
on this sample [2H], we use g = —0.36 and = 55 mT 
vs. 10 mT in the left vs. right dot. Note that dehnes 
the long-time limit of the EDSR frequencies in each dot, 
and hence also tpeak~^peak ia this limit. The knowledge of 
Bnin(r) allows us to attribute the two EDSR resonances 
to the individual dots. We hnd that both, decay times 
and Bnuc(r) generally differ from FP to FP and between 
the two dots, see Fig.j^a). In these measurements, decay 
times range from 3 minutes to almost 8 minutes, much 
longer than earlier hndings in lateral dots isnn]. 

For both types I and III, our measurements show the 
build-up of a strong nuclear polarization which partially 
compensates Rext and Rnm as evidenced by the small val¬ 
ues of Reff, see Fig. [^a). FP I is characterized by a nu¬ 
clear spin polarization which tends to equalize Rgff in the 
two dots: From the small current, we already concluded 
that AR^ ~ 0. From our EDSR measurements, we find 
R^^^ = -(270 ± 5) mT vs. -(320 ± 9) mT in the left vs. 
right dot, corresponding to AR^^ = (—5 ± 10) mT, hence 
ABgff ~ 0. This complete compensation of Bext+Bnm(r) 
is surprising as the EDSR measurement is performed at 
dot positions shifted by ~ 10 nm compared to where 
they are during polarization build-up. It is a first in¬ 
dication that nuclear spin polarization diffuses outside 
the dots. For FP III we find B^^^ = —(340 ± 25) mT vs. 
— (290±20) mT in the left vs. right dot, corresponding to 
ArJjj = (80± 30) mT (with statistical errors from the fit 
procedure). Furthermore, FPs II and IV have similarly 
large leakage currents as FP III, indicating that these 
three FPs have in common a relatively large ABes at 
finite Rext- Unfortunately, we have not been able to reli¬ 
ably capture EDSR current peaks following DNSP traces 
II and IV as these resonances turned out to be too un¬ 
stable. 

In Fig. db), we compare decay curves of the later 
occurring current maximum corresponding to the left 
dot for FP III and three different values of B^^. We 
find equal decay times for relatively small B^'^, namely 
T = 230±20s for R™“ = 0.36 T and 0.65 T, but a longer 
T = 400 ± 60s for 1.48T. This points to a stabiliza¬ 
tion mechanism of the nuclear spin polarization inside 
the dots for longer sweeps. 

To further explore the long decay times observed in 


Fig. 

Fig. 


a) and the stabilization mechanism evidenced in 

b) we continued polarizing with U = 1 mV at a 


fixed R^“ = 0.65 T and measured tpeak and tpgak ^ 


function of the additional polarizing time thoid- As shown 
in Fig. dc) we Hnd that both current maxima occur 
later in proportion to the increase of thoid- Similar life¬ 
time enhancement of B„uc(r) was previously observed in 
bulk GaAs [35] and self-assembled quantum dots [331110]. 
Our long polarization lifetimes of many minutes for large 
thoid and the linear dependence of tpeak on thoid con¬ 
firm a significant polarization of the surroundings of the 
dots. The large and widely spread magnetic field inhomo¬ 
geneity caused by our two nanomagnets further increases 
the lifetime by slowing nuclear spin diffusion via flip-flop 
processes. The latter are restricted by energy conserva¬ 
tion but the nuclear spin state energies depend on the 
local magnetic and electric (and strain) fields, all hav¬ 
ing varying gradients throughout the DQD. In our sys¬ 
tem, the gradient of B^^ is around 0.2-1 mT/nm [35] . 
The ensuing difference in Zeeman energy between closest 
homonuclear atoms is several times larger than their nu¬ 
clear spin dipole-dipole coupling and of the same order as 
the Knight-field gradient caused by the inhomogeneous 
electron wave function and estimated to reduce diffusion 
coefficients by a factor 2-10 m- Since the magnetic field 
gradient extends far beyond the dots [35], it also facili¬ 
tates the polarization of their surroundings. Increasing 
thoid from 0 to 100 minutes causes tpea^ to be delayed 
by more than a factor four while tpe^k only increased 
by 70%. We ascribe this difference to asymmetries in 
geometry and magnetic field in the two dots. 

In summary, we have dynamically polarized nuclear 
spins at a DQD and combined transport spectroscopy 
with electric-dipole-induced spin resonance to study the 
polarization and decay dynamics of nuclear spins. Mea¬ 
suring the leakage current through the DQD in Pauli- 
spin blockade, we And a remarkably complex current 
behavior during magnetic held sweeps. The statistical 
re-occurrence of four patterns in /(Rext) establishes the 
existence of multiple fixed points, one of which is always 
occupied as long as dynamical polarization is maintained. 
Our EDSR measurements reveal long decay times of the 
nuclear spin polarization, its stability being enhanced by 
the strongly inhomogeneous magnetic Held distribution 
generated by two single-domain nanomagnets in an ex¬ 
tended area including the DQD. In addition, the EDSR 
measurements confirm that the individual FPs substan¬ 
tially differ by their polarizations and dynamics of the 
nuclear spins. On one hand, our studies demonstrate 
that the existence of several FPs in dynamical nuclear 
spin polarization complicates the desired control of elec¬ 
tron and nuclear spins in coupled quantum dots. On the 
other hand, our experiments present a salient advance in 
our understanding of the hyperfine induced dynamics in 
nanoelectronic circuits and brings us closer towards the 
desired fine control of the nuclear spins, important for 
quantum information applications. 
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